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Abstract 

The Non-equilibrium Self-consistent Generalized Langevin Equation theory of irreversible relax¬ 
ation [Phys. Rev. E (2010) 82, 061503; ibid. 061504] is applied to the description of the non¬ 
equilibrium processes involved in the spinodal decomposition of suddenly and deeply quenched 
simple liquids. For model liquids with hard-sphere plus attractive (Yukawa or square well) pair 
potential, the theory predicts that the spinodal curve, besides being the threshold of the thermo¬ 
dynamic stability of homogeneous states, is also the borderline between the regions of ergodic and 
non-ergodic homogeneous states. It also predicts that the high-density liquid-glass transition line, 
whose high-temperature limit corresponds to the well-known hard-sphere glass transition, at lower 
temperature intersects the spinodal curve and continues inside the spinodal region as a glass-glass 
transition line. Within the region bounded from below by this low-temperature glass-glass tran¬ 
sition and from above by the spinodal dynamic arrest line we can recognize two distinct domains 
with qualitatively different temperature dependence of various physical properties. We interpret 
these two domains as corresponding to full gas-liquid phase separation conditions and to the for¬ 
mation of physical gels by arrested spinodal decomposition. The resulting theoretical scenario is 
consistent with the corresponding experimental observations in a specific colloidal model system. 

PACS numbers: 23.23,+x, 56.65.Dy 


1 


I. INTRODUCTION. 


Quenching a liquid from supercritical temperatures into the liquid/gas coexistence region 
leads to the separation of the system into the equilibrium coexisting phases [lj, I2I] . This phase 
separation process is a paradigmatic illustration of non-equilibrium irreversible processes, 
whose description must be based on the general fundamental principles of molecular non¬ 


equilibrium statistical thermodynamics [3 


, 4], When the final density and temperature of 


the quench correspond to a state point (n, T ) inside the spinodal region, the separation 


process starts with the amplification of spatial densi 
the two-phase state in thermodynamic equilibrium 


involving mostly colloidal liquids 


11 


y fluctuations, eventually leading to 


10]. However, under some conditions, 


161 ]. this process of phase separation is interrupted 


when the colloidal particles cluster in a percolating network of the denser phase, formin g an 
amorphous sponge-like non-equilibrium bicontinuous structure, typical of physical gels 17]. 

In understanding these experimental observations, computer simulation experiments in 
well-defined model systems have also played an essential role. This includes Brownian dy¬ 
namics (BD) 18] and molecular dynamics (MD) 19] investigations of the density- and 
temperature-dependence of the liquid-gas phase separation kinetics in suddenly-quenched 
Lennard-Jones and square-well 20] liquids, the BD simulation of gel formation in the col¬ 
loidal version of the primitive model of electrolytes (a mixture of equally-sized oppositely 


charged colloids) 12], and MD simulations of 
the two colloidal species of a binary mixture 11 


Digel formation driven by the demixing of 


]. From the theoretical side, on the other 


hand, a fundamental framework is missing that provides an unified description of both, the 
irreversible evolution of the structure and dynamics of a liquid suddenly quenched to an un¬ 
stable homogeneous state, and the possibility that along this irreversible process the system 
encounters the conditions for its structural and dynamical arrest. 

Any theoretical attempt to model this complex phenomenon, however, must be guided 
by the well-established experimental facts and simulation results above, according to which, 
arrested spinodal decomposition is indeed a possible mechanism for gelation in attractive 


colloidal suspensions. Although there is sti 


such as the location of the arrest line 


□ 


1 an ongoing discussion regarding specific details, 


14], one can expect that some of these details 
are actually system-dependent, and that they will eventually be fully understood. In the 
meanwhile, it is important to focus on the simplest experimental model systems, where 
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this phenomenology has been characterized as carefully and as thoroughly as possible. The 
simplest such model system should involve, for example, a perfectly monodisperse suspension 
of spherical particles with attractive interactions that do not require a third species (as 
colloid-polymer mixtures do), so that their interactions are state-independent, without the 
need to translate the real temperature into any form of effective temperature. The solutions 
of the globular protein lysozyme studied by Schurtenberger and collaborators [l3|, [14] are, 
to our knowledge, the closest experimental realization of such an ideal model system. Thus, 
the experimental scenario established by this group may provide an invaluable guide that 
might serve as a reference for the initial theoretical advances in this field. For this reason, 


we now quote the main 
According to Refs. 


ea 


13 


;ures of such experimental scenario. 


14], the experimental equilibrium phase diagram of aqueous 
lysozyme solutions shows at high temperatures a stable gas-crystal equilibrium coexistence, 
and at lower temperatures, a metastable gas-liquid phase separation. Such equilibrium 
scenario is typical of systems of particles interacting through hard-sphere-like repulsions 
plus very short-ranged attractions. The non-equilibrium experiments performed consist of 
the fast temperature quench at fixed overall protein concentration (or volume fraction 0), 
starting with the system in equilibrium at an initial temperature T* above the gas-liquid 
coexistence curve T 9 i (0), cooling the system to a final temperature T lower than 
The region below the coexistence curve “can be separated into three areas differing in their 
kinetic behavior: a region of complete demixing (I), gel formation through an arrested 
spinodal decomposition (II), and an homogeneous attractive glass (III)” [h], Q. More 
specifically, for a shallow quench to a final temperature below the spinodal curve, T < T s (0), 
the classical sequence of spinodal decomposition leading to complete phase separation is 
observed (region I). However, for quenches below a threshold (or “tie-line”) temperature 
T o (0) (< T s (0)), the spinodal domain structure initially coarsens but then completely arrests, 
forming gels via an arrested spinodal decomposition (region II). Microscopically, these gels 
correspond to a coexistence of the dilute fluid with a dense percolated glass phase. At lower 
temperatures, region II is delimited by the boundary with the glass phase (region III), where 
homogeneous attractive glasses can be reached by quenches at sufficiently low T. Many other 
observations, which can be consulted in the original references 12j[ 14], enrich the observed 
experimental scenario with more detail, but the existence of these three kinetically distinct 
regions constitute the most relevant feature calling for a fundamental explanation. 
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Some elements needed to construct the desired unifying theoretical framework may al¬ 
ready be contained in the classical theory of spinodal decomposition js tj. Unfortunately, it 
is not clear how to incorporate in these approaches the main non-equilibrium features of the 
process of dynamic arrest, including its history dependence and aging behavior. Similarly, 
it does not seem obvious how to incorporate the non-stationary evolution of the structure 
and dynamics of a liquid during the process of spinodal decomposition in existing theories of 
glassy behavior [21]. For example, although many of the predictions of conventional mode 
coupling theory (MCT) 22-1251] have found beautiful experimental verification, this theory 
is unable to describe irreversible non-stationary processes. The reason for this is that, in its 
present form, MCT is in reality a theory of the dynamics of liquids in their thermodynamic 


equilibrium states. To overcome this fundamental limitation, in 2000 Latz 


26] proposed 


a formal extension of MCT to situations far away from equilibrium which, however, has 
not yet found a specific quantitative application. In this context, let us also refer to the 
mean-field theory of the aging dynamics of glassy spin systems developed by Cugliandolo 
and Kurchan 27]. This theory has also made relevant detailed predictions that have been 
verified in experiments and simulations. Unfortunately, the models involved lack a geometric 
structure and hence cannot describe the spatial evolution of real colloidal glass formers. 

MCT’s limitation above is shared by the self-consistent generalized Langevin equation 


(SCGLE) theory of dynamic arrest 


3-30|, 


a theory that predicts, in a manner completely 


analogous to MCT, the existence of dynamically arrested states. In contrast to MCT, how¬ 
ever, the SCGLE theory has been extended to describe the irreversible non-equilibrium 
evolution of glass-forming liquids, thus resulting in the “non-equilibrium self-consistent gen¬ 
eralized Langevin equation” (NE-SCGLE) theory [31M33I]. This extended theory provides a 
conceptually simple picture of the crossover from ergodic equilibration to non-equilibrium 


aging 


34|[35|], offering a perspective of the glass transition in which the “waiting”-time t be¬ 


comes a relevant active variable and in which indeed the measured properties of the system 
may depend on the protocol of preparation and on the duration of the actual measurements. 
Thus, this non-equilibrium theory may provide the general framework in which to build the 
missing unified description of the processes of arrested spinodal decomposition. 

The main purpose of the present paper is to take the first steps in exploring this possibility. 
The first of these steps is to discuss the manner in which the NE-SCGLE theory may be 
adapted to address the description of arrested spinodal decomposition. The second refers to 
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obtaining and presenting the concrete results of its application to a specific model system 
that permits at least qualitative contact with the observed behavior of an experimental model 
system such as that described above. The physical interpretation of our theoretical results, 
however, is not as straightforward as, for example, calculating the equilibrium properties of 
a system using century-old and well-established statistical thermodynamic methods (such 


as calculating the partition function Q). Instead, because of the non-equilibrium and non¬ 
linear nature of the phenomenon described, the interpretation of the specific predictions is 
much more subtle, involving features that might seem at first sight strongly counterintuitive. 
Thus, the third and most relevant step, is to properly interpret these results. 


Regarding the first of these steps, let us mention that the mos 


NE-SCGLE theory of irreversible processes in colloidal liquids 31 


general version of the 


33] consists of two fun¬ 


damental time-evolution equations, one for the mean value n(r,f), and another for the 


covariance a(r,r / ;t) = Sn(r, t)5n(r', t) of the fluctuations 5n(r,t) = n(r,t) — n(r,f) of the 
local concentration profile n(r, t ) of a colloidal liquid. Thus, in the following section and in 
the appendix we review this general description, to discuss the form adopted by these equa¬ 
tions in the absence of external fields, and the possibility that they allow the description of a 
fingerprint of spinodal decomposition, namely, the early stage of the amplification of spatial 
heterogeneities. The resulting equations are thus restricted to allow the discussion of the 
simplest explicit protocol of thermal and mechanical preparation, namely, an instantaneous 
isochoric quench, to describe the spontaneous relaxation of the system toward its equilibrium 
state or toward predictable arrested states. To simplify the reading, Section |TT] is essentially 
a summary of the specific NE-SCGLE equations that we shall actually solve, while the ap¬ 
pendix contains a description of the main derivations and approximations leading to these 
specific equations. 

In Sect. m we apply the resulting approximate theory to the discussion of the possi¬ 
bility that the process of spinodal decomposition may be interrupted by the emergence of 
dynamic arrest conditions. For the sake of concreteness, there we shall have in mind a very 
specific model system, namely, the “hard-sphere plus attractive Yukawa” (HSAY) potential, 
subjected to an instantaneous quench at t — 0 from an initial temperature Tj to a final 
temperature Tf = T (< Tj), with the volume fraction </> held fixed. The main focus of this 
application is the determination of the dynamic arrest diagram in the region of unstable ho¬ 
mogeneous states. According to the predicted scenario, the spinodal curve turns out to be 
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the borderline between the regions of ergodic and non-ergodic homogeneous states, and the 
high-density liquid-glass transition line, whose high-temperature limit corresponds to the 
well-known hard-sphere glass transition, actually intersects the spinodal curve at lower tem¬ 
peratures and densities, and continues inside the spinodal region as a glass-glass transition 
line. 

Some elements of this predicted scenario, however, may seem strongly counterintuitive, 
and not consistent with experimental observations. For example, this scenario predicts that 
any quench inside the spinodal region will meet conditions for dynamic arrest, while we know 
that at least for shallow quenches the system will definitely phase-separate. A satisfactory 
understanding of these features requires a detailed discussion of a number of issues whose 
subtlety derives from the lack of familiarity with the solutions of the full NE-SCGLE non¬ 
linear equations. Thus, in Sect. |IV]we discuss the physical interpretation of this theoretical 
scenario by analyzing in more detail the predicted long-time non-equilibrium arrested struc¬ 
ture factor S a (k; 0, T), particularly the length scale £ a (0, T) = 2ti/ k max associated with the 
position k max of its small-/;; peak. This length represents the size of the spatial heterogeneities 
typical of the early stage of spinodal decomposition described by the Cahn-Hilliard-Cook 
(CHC) classical linear theory of spinodal decomposition [ 5 , 6], shown to be a particular 
limit of the NE-SCGLE theory. This analysis results in the recognition that indeed, shallow 
quenches lead to full phase separation whereas deeper quenches lead to the formation of gels. 
Thus, in Sect. we hnally establish a more detailed connection between our theoretical 
scenario and the experimental observations in the lysozyme real model system, and in Sect. 
EO we summarize the main results of this work. 


II. REVIEW OF THE NE-SCGLE THEORY. 

Let us start by considering a Brownian liquid formed by N particles that execute Brown¬ 
ian motion (characterized by a short-time self-diffusion coefficient Dq) in a volume V while 
interacting between them through a generic pair potential u(r) = UHs ( r ) + ua(j') that is the 
sum of a hard-sphere (HS) term uhs { t ) with HS diameter a, plus an attractive tail ua{t). 
Let us imagine that we subject this system to a prescribed protocol of thermal and/or me¬ 
chanical treatment, described by the chosen temporal variation of the reservoir temperature 
T(t), of the total volume V(t), and of the applied external fields 0(r, t). In principle, the 
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ultimate goal of the NE-SCGLE theory is to predict the response of our system to each pos¬ 
sible protocol of thermal and mechanical manipulation. For simplicity let us assume that 
the thermal diffusivity of the system is sufficiently large that any temperature gradient is 
dissipated almost instantly (compared with the relaxation of chemical potential gradients), 
so that the temperature held inside the system can be considered uniform, and instantly 
equal to the temperature of the reservoir, T(r,t) = T(t). 

Let us now imagine that we use these time-dependent fields and constraints to drive the 
system to a prescribed (but arbitrary) initial state characterized by a given mean concen¬ 
tration profile n°(r) and covariance a°(fc; r), and we then fix the held, the temperature, and 
the volume afterward, i/)(r,t) = ^(r), T(t) = T, and V(t) — V for t > 0. The NE-SCGLE 
theory will then allow us to predict the most fundamental and primary information that 
characterizes the intrinsic properties of the system, namely, its spontaneous behavior during 
the irreversible relaxation towards its new expected equilibrium state. This is described in 
terms of the time-evolution of the non-equilibrium mean concentration prohle n( r, t) and 
covariance a(k] r ,t) for t > 0 (see Eqs. (P and (J2J) of the appendix). 

In Ref. [34j we started the systematic discussion of this topic in the context of a model 
system involving only purely repulsive interactions. For this, we restricted ourselves fur¬ 
ther to the simplest explicit protocol of thermal and mechanical preparation of our system, 
namely, its instantaneous cooling or heating from an initial temperature T t to a final tem¬ 
perature Tf, so that T{t) = Tj for t < 0 and T(t) — Tf for t > 0, under isochoric conditions, 
V{t) = V (so that n{t) = n) at all times, and in the absence of applied external fields, 
-0(r, t) = 0. As explained in Ref. 34J], the NE-SCGLE theory will then describe the spon¬ 
taneous response of the system in terms only of the non-stationary but uniform covariance 
a[k\ t), now written as the non-equilibrium static structure factor S(k] t ) = a[k\ t)/n, whose 
time-evolution equation for t > 0 reads 

dS{k\t) 


dt 


= —2k 2 D^b(t)n£f(k) [S(k\t ) — 1 /nSj(k)\ 


( 2 . 1 ) 


In this equation £f(k) = Safari, Tf) is the Fourier transform of the thermodynamic 
functional derivative £[r,r';n,T\ = [(f/3/i[r; n, T\/5n(r')], evaluated at the uniform den¬ 
sity and temperature profiles n(r) = n and T( r) = T, in which case it can be written 
as £[r, r'; n, T] = £ h {|r — r'|; n, T ) = <5(r — r ')/n — c(|r — r'|; n, T ) or, in Fourier space, as 
£h(k\ n, T) = 1/n — c{k\ n,T), where c(r; n, T) is the so-called direct correlation function. 
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The non-stationary, non-equilibrium mobility function b(t) entering in Eq. (12.ip above, 
is the most important kinetic property. Within the NE-SCGLE theory this property is 
determined by (see Refs. 32l, l33f) 


b(t) = [1 + / drAC(r-,t)] , 

Jo 

with the t-evolving, r-dependent friction coefficient Agiven approximately by 


( 2 . 2 ) 


A C*(r;t) = 


D, 


dk k 


S(k-t) - 1 
S(k-t) 


(2.3) 


247r 3 n J 

xF(k,T]t)F s (k,T]t) 

in terms of S(k; t) itself and of the collective and self non-equilibrium intermediate scattering 
functions F(k, r; t ) and Fs(k, t; t ), whose memory-function equations are written, in Laplace 
space, as 

. S(k-t) 


and 


F s (k,z]t) = 


z + 

k 2 D°S~ 1 (k;t ) ’ 

1+A(fc) A 


1 

z + 

k 2 D° 


(2.4) 


(2.5) 


with F(k,z]t), F s (k,z]t ), and A (*(z]t), being the corresponding Laplace transforms (LT) 


and with A (k) being a phenomenological “interpolating function” 30 1 ], given by 

A (k) = 1/[1 + {k/k c 


vc) 2 ], 


( 2 . 6 ) 


where the phenomenological cut-off wave-vector k c depends on the system considered. Since 
in this paper we are only interested in semi-quantitative trends, here we adopt the value 
k c = 1.305(27t/cj) = 8.2/cr, determined in a calibration procedure involving simulation data 


of the hard-spheres system 


4l|. 


III. DYNAMIC ARREST DIAGRAM OF A SIMPLE MODEL SYSTEM. 


In this section we apply the previous general theory to the discussion of the possibility 
that the process of spinodal decomposition may be interrupted by the emergence of dynamic 
arrest conditions. For this, let us consider a generic liquid formed by N particles in a volume 


















V interacting through a simple pair potential u{r) = Uns ( r ) + u a(t), which is the sum of 
a hard-sphere term uhs(j) plus an attractive tail ua(t). Although all the results that we 
shall discuss here are independent of the specific functional form of ua(t), for the sake of 
concreteness, in what follows we shall have in mind a very specific model system, namely, 
the “hard-sphere plus attractive Yukawa” (HSAY) potential, defined as 


u(r) 


oo, 

r < o 

,exp[— z(r/<7— 1)] 

r > a, 

(r/cr) 


(3.1) 


where a is the hard sphere diameter, e the depth of the attractive Yukawa well at contact, 
and z^ 1 its decay length (in units of a). For given er, e, and z, the state space of this system 
is spanned by two macroscopic variables, namely, the number density n = N/V and the 
temperature T, which we express in dimensionless form, as n* = no 3 and T* = k B T/e. In 
practice, however, from now on we shall use o as the unit of length, and e/k B as the unit 
of temperature, so that n* = n and T* = T; most frequently, however, we shall refer to the 
hard-sphere volume fraction cj) = nn/Q. 


A. Equilibrium thermodynamic properties. 


The most fundamental thermodynamic property of this system, that we need to specify 
in order to apply the NE-SCGLE theory described above, is the thermodynamic functional 
derivative £[r,r';n,T\ = [<J/fyt[r; n, T\/Sn(r')\. Evaluated at the uniform density and tem¬ 
perature profiles n(r) = n and T(r) = T, we have that £[r, r';n, T] = Sh {|r — r'|; n,T) = 
5(r — r')/n — c(|r — r'|; n,T) or, in Fourier space, as £h(k]n,T) = 1/n — c(k]n,T), with 
c(r;n,T) being the direct correlation function. In the present paper we shall rely on the 
approximate prescription proposed by Sharma and Sharma 42] to determine this thermo¬ 
dynamic property, which for the HSAY potential above reads 


c(k\ <t>, T) = c H3 (k\ <t>) - fiu A (k), 


(3.2) 


in which is the FT of the direct correlation function of the hard-sphere 

liquid, approximated by the Percus-Yevick approximation with Verlet-Weis correction 
(PYVW) 43, M\ and (3iiA{k) is the FT of the attractive potential (in units of the 
thermal energy k B T = 1//3), i.e., /3uA(k) = An drr 2 [sm(kr)/(kr)]/3uA(r). Thus, 
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FIG. 1: Binodal and spinodal lines of the gas-liquid coexistence of the HSAY model fluid in 
Eq. (13.11) with z = 2, calculated with the Sharma-Sharma approximation for £f l (k',n,T). The 
corresponding critical point (black star) is located at (0, T) = (0.13,0.85). Also shown is the 
freezing line, calculated from the Hansen-Verlet condition that the maximum of S^ eq '(k] 0, T) = 
l/n£h(k] 0, T) reaches 2.85. The corresponding triple point (soft red circle) is located at (0, T) = 
(0.33,0.65). Finally, the dark blue solid line is the liquid-glass dynamic arrest line, calculated from 
the solution 7 eq (T, fa) of the equilibrium “bifurcation equation”, Eq. (13.31) . The dark circle denotes 
the intersection (fa, fa) = (0.33,0.45) of this arrest line with the spinodal curve. 

for the HSAY potential, (3u a (k) = —(4n/T*) f™ drr 2 [sin(kr)/(kr)] exp[— z(r — l)]/r = 
— (47 t/T*) [( kcosk + zsin k)/k(k 2 + z 2 )]. 

For fixed z, the equilibrium phase diagram of the HSAY model system in the state space 
(cf), T ) contains the gas, liquid and (crystalline) solid phases. The previous approximation 
for £h(k',n,T) allows us to sketch the most prominent features of the fluid phases. For 
example, the gas-liquid transition involves a coexistence region, with its associated binodal 
and spinodal lines. The former can be determined by integration of fadp/drifa = Efak = 
0; (j), T ) to get the equation of state, along with Maxwell’s construction. The latter can be 
determined from the condition for thermodynamic instability of uniform fluid states, namely, 
£h(k = 0; n , T) < 0. The resulting binodal and spinodal lines are illustrated in Fig. |Tj for 
the HSAY system with z = 2. On the other hand, the equilibrium static structure factor of 
the fluid state is determined by the equilibrium (Ornstein-Zernike) condition S^ eq \k\ (f>, T) = 
l/n£h(k] (f>, T ). From S^ eq \k; 0, T ) the freezing line, which is another boundary of stability of 
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the uniform 
condition 


iquid state, can be readily sketched using the phenomenological Hansen-Verlet 
45j that the height Smax{4> 1 T) °f the main peak of the equilibrium static structure 


factor S( eq \k;(/),T) reaches the value Smal{4>,T) = 2.85. Besides the spinodal, binodal and 
freezing lines, Fig. |T| also exhibits the corresponding gas-liquid critical point, located at 
(0, T) = (0.13,0.85), and the triple point at (0, T) = (0.33,0.65). These equilibrium lines 
serve as a reference to the introduction, in the following subsection, of the liquid-glass 
dynamic arrest line, also obtained from S( eq \k-,(j),T) = l/n£h(k;(p,T). 


B. Dynamic arrest diagram from the equilibrium SCGLE theory. 

The first relevant task in the description of dynamic arrest phenomena in a given sys¬ 
tem, such as that in our example, is the determination of the dynamic arrest diagram, i.e., 
the region containing the state points (0, T) where the system will be able to reach ther¬ 
modynamic equilibrium (ergodic region) and the region where the system, prevented from 
crystallizing, will be trapped in a dynamically arrested state (non-ergodic region). Under 
normal circumstances this can be approached solving the so-called bifurcation equations 


of mode coupling theory 


22 ], or the corresponding equations of the equilibrium version of 


the SCGLE theory 30]- The latter consists of the following equation for the equilibrium 
localization length squared 7 ^ eiJ q0, T), 


dkk 4 


[S^\k,4>,T) - l ] 2 A 2 (fe) 


(3.3) 


y( e( ?) 67 r 2 n J 0 [A(fc)S'( e 'b(/c; 0, T) + k 2 ^ eq l\ [A(fc) + k 2 ^^] 

The only input of this equation is the equilibrium static structure factor S ( ' eq \k; 0, T) = 
1 /n£h(k;4>,T) at a given state point (0, T). If the solution is 7 ^(0,T) = 00 , we conclude 
that the state point ( 0 , T ) lies in the ergodic region, whereas if j( eq \(j), T ) is finite, the point 
(0, T) lies in the non-ergodic region of state space. I 11 this manner we can determine the 
dynamic arrest transition line. 

This procedure has been performed in the past for several model systems, for which the 
equilibrium static structure factor S^ eq \k\ 0, T), given in terms of the equilibrium condition 
Ste>(M,T) = l/n£ h (fc;0,T), is well-defined in the entire state space (0, T ). In the present 
case we have carried out this exercise for our HSAY model, and the result is the liquid- 
glass dynamic arrest transition line (0 C , T c ) represented in Fig. [I] by the dark (blue) solid 
line, with the region to the left and above this curve corresponding to ergodic states. This 
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dynamic arrest line extends from its high temperature limit, corresponding to the hard- 
sphere glass transition at (0, T) = (0.582, oo), down to its intersection (0b, Tfi) = (0.33, 0.45) 
with the spinodal curve, indicated by the empty circle in the figure. Unfortunately, for 
systems with thermodynamically unstable (spinodal) regions, like our HSAY model, we can 
only apply this method outside such unstable region, since for the state points (0, T) inside 
the spinodal no equilibrium static structure factor S^(k',<j>,T) exists that corresponds to 
spatially uniform states. It is at this point that the power of the non-equilibrium version of 
the SCGLE theory is manifested, since this general theory does provide a proper manner to 
overcome this limitation of the equilibrium theory thus offering an intriguing and unexpected 
qualitative scenario of the interplay between dynamic arrest and gas-liquid phase separation. 
As we shall see below, the alternative method is NOT based on the use of the equilibrium 
bifurcation equation, Eq. (13.3[) , but on its non-equilibrium generalization, Eq. (13.7ft below. 

C. Detecting dynamic arrest transitions using the ATE-SCGLE theory. 

Let us now use the NE-SCGLE theory to detect dynamic arrest transitions. We start 
by noticing that from the very structure of Eq. (ED one can infer the existence of two 
fundamentally different kinds of stationary solutions, representing two fundamentally dif¬ 
ferent kinds of stationary states of matter. The first corresponds to the condition in which 
stationarity is attained because the factor [S(k] t) — l/nSf(k)\ on the right side of Eq. 
(12.ip vanishes, i.e., because S(k]t ) is able to reach its thermodynamic equilibrium value 
S^'farijTf) = 1 /n£h(k',n,Tf), while the mobility b(t ) attains a finite positive long-time 
limit by. These stationary solutions correspond, of course, to the ordinary thermodynamic 
equilibrium states, in which the difference [S(k; t)—l/nSf{k)\ decays (according to Ec. (12.11) . 
and in a gross approximation) exponentially fast, [S(k]t) — 1 /n£f(k)\ oc exp[— t/t eq ], with 
the equilibration time t eq estimated as t eq (n,Tf ) oc [2k^ nax D°n£f(k max )b e y] \ 

The second class of stationary solutions of Eq. (12. ip emerges from the possibility that the 
long-time asymptotic limit of the kinetic factor b(t) vanishes, so that dS(k ; t)/dt vanishes at 
long times without requiring the equilibrium condition [S(k;t) — 1 /n£ h {k-,n,Tf)\ = 0 to be 
fulfilled. This mathematical possibility has profound and general implications. For exam¬ 
ple, under these conditions S(kfi) will now approach a distinct non-equilibrium stationary 
limit, denoted by S a (k ), which is definitely different from the expected equilibrium value 
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S (eq ^(k: n,Tf) = l/n£ h (k;n,Tf). Furthermore, the difference [S(k;t) — S a (k)} is predicted 
to decay to zero not exponentially, but in a much slower fashion, namely, as t~ x , with an 
exponent x near unity [331, |34j. Contrary to ordinary equilibrium states, whose properties 
are determined solely by the fundamental condition of maximum entropy, the properties of 
these stationary but intrinsically non-equilibrium states, such as S a (k), may depend on the 
preparation protocol (in our example of the instantaneous isochoric quench, on T) and Tf). 

For a given instantaneous isochoric quench from an initial temperature 7) to a lower final 
temperature Tf, one can also predict if the system will equilibrate or will be trapped in 
a non-equilibrium arrested state by analyzing the stationary solutions of Eq. (12.11) . This 
analysis simplifies greatly with the change of variables from the actual evolution time t to 
the so-called “material” time u(t) defined by 34] 


u(t) = / b(t')dt'. 


(3.4) 


This allows us to write, for example, the actual solution S(k\t ) of Eq. (12.11) as S(k',t) = 
S*(k;u(t)), with the function S*(k',u) being the solution of 

dS*(k;u ) 


du 


= -2 k 2 D"n£ f (k) [S*(k]u ) - 1 /fiS f (k)\ 


(3.5) 


i.e., 


S(t,t) = S*(k; u(t)) = [nSfik)}- 1 + {Si(k) - [nE^k)}- 1 } e - 2 * 2 A>»*/(*)«(*), ( 3 . 6 ) 


where Si(k) = S(k]t = 0 ) = S*(k;u = 0 ) represents the (arbitrary) initial condition. 

This expression interpolates S*(k\u ) between its initial value S^k) and the value of the 
thermodynamic property \n£f(k)]~ 1 . Thus, if nothing impedes the system from reaching 
equilibrium at the state point (n, Tf), the non-stationary static structure factor S(k]t ) will 
eventually attain its equilibrium value S^ eq \k;n,Tf) = \n£f(k)]~ 1 . This, however, is only 
one of the two possibilities described above. In order to determine which of them will 
govern the course of the spontaneous response of our system, we actually input S*(k;u ) in 
the “bifurcation” equation for the square localization length y(t) at evolution time t (i.e., 
the long-r asymptotic value of the mean squared displacement). Denoting y(t) = 7 *{u{t)) 
simply as 7 *(w), such an equation reads 


_j_ = _j_ r ,,,4_ is'(t,u)-ifx 2 (t,u) _ 

7 *(u) 67 r 2 n J 0 [X(k]u)S*(k]u) + k 2 'y*(u)][X(k]u) + k 2 'y*(u)] 


(3.7) 
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As explained in Ref. 34], if we find that 7 *(w) = oo for 0 < u < oo, we conclude that the 
system will be able to equilibrate after this quench, and hence, that the point ( n,Tf ) lies in 
the ergodic region. If, instead, a finite value u a of the parameter u exists, such that 7 *{u) 
remains infinite only within a finite interval 0 < u < u a , the system will no longer equilibrate, 
but will become kinetically arrested, with the non-equilibrium asymptotic structure factor 
S a {k) given by 


S a (k ) = S*{k;u a ) = [nSfik)]- 1 + {Si(k) - [n^)]” 1 } e ~^on£ f {k)u a _ ( 3 . 8 ) 


The parameter u a and the corresponding long-time asymptotic square localization length 
7* = 7 *( u a) are dynamic order parameters in the sense that if u a = 7* = cx), the system will 
reach equilibrium, but if u a and 7* are finite, the system will be dynamically arrested. Both 
of them, as well as the non-equilibrium static structure factor S a (k), will depend on the 
protocol of the quench (in our case, on </>, Tj, and the final temperature Tf, which from now 
on will be denoted simply as T ). Thus, if one determines the functions u a = u a {(ft, Tj,T), 
7* = 7*(0, Tj,T), and S a (k) = S a (k] <f), T), T), one can in principle draw the dynamic arrest 
diagram. 

We have implemented this protocol and here we present the results of its application to 
the concrete illustrative case involving the HSAY model. Let us first consider isochores with 
(j) > (j>b , which lie to the right of the intersection indicated by the dark circle in Fig. [D A 
representative isochoric quench of this class is schematically indicated in the inset of Fig. 
Ha) by the downward arrow along the isochore (j) = 0.4. In the main panel of this figure 
we plot the resulting m“ 1 (T) as a function of T (solid line), together with the corresponding 
value of 7a” 1 (T) (dashed line). 

The first feature to notice is that indeed, for quenches to final temperature T above the 
critical temperature T c = 0.702, u a (T ) is always infinite (so that u^iT) = 0, and hence, do 
not appear in the logarithmic window of the figure), whereas for T <T C we detect a finite 
value of u a (T). This value, however, diverges as T approaches T c from below, and decreases 
monotonically as the final temperature T decreases. According to the results in Fig. [ 2 ]( a )> 
for T <T C also 7*(T) attains finite values, thus reaffirming the existence of an arrested phase 
in this temperature regime. In fact, these results for 7*(T) reveal a second relevant feature, 
namely, that this dynamic order parameter changes discontinuously at T = T c , from a value 
7* _1 (T“) = 41 to a value 7* -1 (T c + ) = 0. We notice that this jump in the value of 7*(T) 
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FIG. 2: (a) Dependence of u a (T ) (solid line) and of 7 *(T) (dashed line) on the final temperature T of 
an instantaneous quench with initial temperature T = 1.0 and fixed volume fraction <f> = 0.4. The 
inset reproduces essentially Fig. [T] with the arrow indicating this illustrative isochoric quench, (b) 
u a (T ) for the previous quench and for other two quenches that start at different initial temperatures, 
namely, Tj = 1.3 and with Tj = 0.8. Notice that the three cases yield the same result for the 
transition temperature, T c (0 = 0.4) = 0.7. 

corresponds to the jump in the solution 7 eg (T) of the equilibrium bifurcation equation, Eq. 
(13.31) . and it occurs at the same critical temperature T c = 0.702, thus confirming that both 
methods refer to the same dynamic arrest transition, and that this is a “type B” dynamic 
arrest transition in the language of mode coupling theory. To complete this identification, 
we repeated this procedure at other volume fractions to determine the ^-dependence of the 
dynamic arrest transition temperature T c (0), thus finding that indeed, outside the spinodal 
curve (i.e., for volume fractions 0 larger than the volume fraction of the intersection point 
(0b, T b ) = (0.33,0.45)) this method leads to exactly the same liquid-glass dynamic arrest 
line previously obtained using Eq. (13.31) . 

Let us notice, however, that the curves describing the dependence of u a (T ) on the final 
temperature T do depend on the initial temperature Tj of the quench. However, the resulting 
transition temperature, T c (0 = 0.4) = 0.702, is independent of Tj. This is illustrated in Fig. 
|2](b), which compares the function u a (T) obtained with Tj = 1.0 with the results for u a (T ) 
obtained with Tj = 1.3 and Tj = 0.8. This independence of T c (0) on the initial temperature 
Tj is also a necessary condition for the present method to be a reliable approach to the 
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determination of the liquid-glass dynamic arrest transition. In Fig. [21(b) we have also 
included the results for y* _ 1 (T) corresponding to the quench processes with different initial 
temperatures, to show that 7 *(T) is much less sensitive to the initial temperature of the 
quench. Let us also notice that for 0 > 0 6 , the results for 7 *(T) and u a {T ) as a function 
of the final temperature T, do not exhibit any discontinuous behavior when crossing the 
spinodal line, i.e., their dependence on T does not detect the spinodal. Thus, for 0 > 0b the 
spinodal condition £0(fc = 0; 0, T) = 0 does not seem to have any significant effect on the 
dynamics of the system. As discussed below, this rather astonishing result is a consequence 
of the structure of the expression for S a (k ) in Eq. (I3.8jb and is in dramatic contrast with 
what happens in the complementary regime, i.e., for volume fractions 0 < 0b, as we now 
explain. 

D. Dynamic arrest inside the spinodal region. 

The main advantage of the NE-SCGLE methodology just explained, for determining the 
dynamic arrest transition line, is that its use can be extended to the interior of the spinodal 
region. This subsection continues the description of the use of this methodology, but now 
to unveil the structure of the dynamic arrest diagram of our model in the region inside the 
spinodal line. Thus, let us repeat the calculations illustrated in Fig. El but now for isochores 
corresponding to volume fractions 0 smaller than the volume fraction 0b = 0.33 of the 
intersection between the liquid-glass dynamic arrest transition and the spinodal curve. The 
corresponding illustrative example is contained in Fig. [3l which plots the functions u a (T) 
and 7 *(T) as a function of T for the HSAY model system subjected to an instantaneous 
quench from an initial temperature T t = 1.0 to a final temperature Tf < 1.0 at fixed 0 = 0.2 
(downward arrow in the inset of the figure). 

In contrast with the results for the isochore 0 = 0.4, in Fig. [21 we notice that in the present 
case the functions u a (T ) and 7 *(T) remain infinite for all temperatures above a critical 
temperature T s = 0.76, and that this singular temperature coincides precisely with the 
temperature of the spinodal curve for that isochore. The results in Fig. [3] also indicate that 
both parameters, u a (T ) and 7 *(T), are finite for T < T s , and both diverge as T approaches 
T s from below. This reveals a rather dramatic and unexpected conclusion, namely, that 
the spinodal curve, besides being the threshold of thermodynamic instability, now turns out 
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FIG. 3: Dependence of u a (T) (solid line) and of 7 *(T) (dashed line) on the final temperature T of 
an instantaneous quench with initial temperature T) = 1.0 and fixed volume fraction 0 = 0.2. The 
inset reproduces essentially Fig. [T] with the arrow indicating this illustrative isochoric quench. 
Notice that y a and u a are discontinuous at a temperature T c , which implies the existence of a 
type-B glass-glass transition, whereas as T approaches to T s 7 a and u a diverge, thus implying a 
type-A dynamic arrest transition at the spinodal. 

to be also the threshold of non-ergodicity. Furthermore, the fact that 7 *(T) diverges as T 
approaches T s from below determines that this transition from ergodic to non-ergodic states 
occurs in a continuous fashion, i.e., that it is a “type A” dynamic arrest transition in MCT 
language. 

Examining again the same results in Fig. [31 but now at temperatures well below the 
spinodal curve, we see that the parameters u a (T ) and 7 *(T) exhibit a discontinuity at a lower 
temperature T c (</>). This discontinuity reveals the existence of still a second dynamic arrest 
transition, now corresponding to a glass-glass “type B” transition, in which the dynamic 
order parameter 7 *(T) changes discontinuously by about one order of magnitude, from a 
value 7 a(T c _ ) = 36 to another finite value 7 q(T c + ) = 4. Notice that the value y*(T _ ) is 
similar to the corresponding value at the isochore (j) = 0.4. Repeating these calculations at 
other isochores, we determine both transition temperatures, T c (</>) and T s (0), as a function 
of volume fraction. The corresponding dynamic arrest transition lines are presented in Fig. 
[4] This figure reveals a rather unexpected global scenario, in which (a) this low-T (type 
B) glass-glass transition (dark blue dashed curve) turns out to be just the continuation 
to the interior of the spinodal region, of the liquid-glass transition previously determined 
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FIG. 4: Dynamic arrest diagram exhibiting the low-T (type B) glass-glass transition (dark blue 
dashed curve), which continues the liquid-glass transition (dark blue solid curve) to the interior of 
the spinodal region. The (red) dotted line that superimposes on the spinodal curve is a continuous 
(type A) ergodic-nonergodic transition. These two transitions merge at the bifurcation point 
(<j>b, Tb). The soft blue dashed curve to the right of the bifurcation point (06, Tb) is the dynamically 
irrelevant part of the spinodal curve (see Subsection IIV E|) . 


outside the spinodal (dark blue solid curve), (b) the spinodal line itself is a continuous (type 
A) ergodic-nonergodic transition, (c) these two transitions merge at the bifurcation point 
( 06 , Tb) to continue outside the spinodal as the liquid-glass transition line that terminates at 
the hard-sphere glass transition point (0, T) = (0.582, oo), and (d) for 0 > 0 b the spinodal 
line does not have any apparent dynamic significance. 

At first glance some of these predictions might appear counterintuitive. Let us remind, 
however, that they are based only on the analysis of the dependence of the dynamic or¬ 
der parameter 7 *(0, T) on the final temperature T of the instantaneous isochoric quenches 
analyzed. This order parameter is a functional of the long-time asymptotic value S a (k ) of 
the non-stationary static structure factor S(k]t). Thus, the analysis of the dependence of 
S a (k] 0, T) on 0 and T should provide additional pieces of structural information, while the 
waiting-time dependence of S(k\t ) (and of all the other structural and dynamic properties) 
will provide valuable pieces of kinetic information. When these pieces are assembled together 
in a structural and kinetic jigsaw puzzle, one should expect that a sound and convincing 
scenario will emerge. For this, however, we must first identify and describe the most rele¬ 
vant of these pieces. The dynamic arrest diagram in Fig. [4] constitutes the initial and most 
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basic piece for this discussion, but in the rest of the paper we shall identify others, equally 
fundamental. 

IV. PHYSICAL INTERPRETATION OF THE DYNAMIC ARREST DIAGRAM. 

The scenario described by the dynamic arrest diagram in Fig. [4] is highly provoking, but 
also at first sight difficult to reconcile with experience. First, we all know that simple fluids, 
when quenched to the inside of the spinodal region, are expected to fully phase-separate 
into the gas phase and its corresponding condensed (liquid or crystal) phase. On the other 
hand, as explained in the introduction, there are numerous experimental and numerical 
evidences that for sufficiently deep quenches, the process of spinodal decomposition may 
be interrupted by its interference with the dynamic arrest of the condensed phase. Thus, 
the dynamic arrest diagram in Fig. [4] may contain the fundamental explanation of this 
experimental double scenario of full vs. arrested spinodal decomposition. 

To test this expectation, however, we must first understand what exactly is the detailed 
physical meaning of these NE-SCGLE predictions and what are their most relevant and 
verihable physical consequences. In this section we discuss several pieces of information 
that derive from the solution of the NE-SCGLE equations, and which will dissipate some 
of the most subtle puzzles of the interpretation of this dynamic arrest diagram. With these 
details taken into account, we shall see that the resulting theoretical scenario is actually 
quite consistent with the main features of the experimental scenario of arrested spinodal 
decomposition in lysozyme protein solutions, thus establishing a concrete link of our theory 
with physical reality. 

A. S(k;t), S^Xk-,(j),T), and S a (fc) below the spinodal line. 

In this subsection, for example, we discuss the mechanism that allows the existence 
of a well-behaved non-equilibrium S(k’,t) for quenches to state points (</>, T) inside the 
spinodal region, where the uniform fluid state has become thermodynamically unstable and 
the equilibrium structure factor S^ eq \k\(j),T) does not exist. According to Eqs. (13.61) and 
(13.8(1 . the calculation of S[k ; t) and S a (k) requires the value of the thermodynamic property 
£h(k ; </>, T ) for k > 0. As we know, the condition £ h {k = 0; 0, T) > 0 is a condition for the 
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stability of uniform thermodynamic states, a condition held for all temperatures T above 
the spinodal temperature T s (<ft). Thus, the spinodal condition £h(k = 0; 0, T s ) = 0 defines 
the threshold of thermodynamic instability of uniform states, so that for all states with 
temperatures below T s (cft) we must have that £h(k ; (ft, T) < 0 not only at k — 0, but also at 
least within a finite interval 0 < k < k 0 . 

This means that the equilibrium static structure factor (k; (ft,T) = l/n£h(k; (ft,T) 
will attain unphysical negative values in this interval and will exhibit a singularity at k — 
k(). This unphysical behavior, which is a manifestation of the non-existence of spatially 
uniform equilibrium states inside the spinodal region, is illustrated in Fig. 0 which plots 
[n£h(k] (ft, T)] -1 (dashed curves) at the state points ((ft, T) = (0.2, 0.2) and (0.4, 0.2) below 
the spinodal curve. One might then expect that the non-equilibrium static structure factor 
S(k\t ) and its asymptotic value S a (k), which are formally an exponential interpolation 
between Si(k) and \n£f(k] (ft, T)] -1 (see Eqs. (13.6p and (13.81) 1. will inherit these singular 
features of the function [n£f(k] (ft, T)] _1 . 

In reality, this is not the case, since the appearance of \n£f(k)\ in the exponent of the 
interpolating function, allows us to show, for example, that near k = ko we have that S(k; t ) 
and S a (k) can be approximated, respectively, by S(k;t) = S*(k-,u(t )) ~ Si(k) + 2 k 2 D 0 u(t) 
and S a (k ) ~ Si(k) + 2k 2 D 0 u a , which are always positive definite. This is illustrated in Fig. [5] 
with the results for S a (k', (ft, T ) (solid lines) for the quenches along the isochores (ft = 0.2 and 
0.4, from the same initial temperature T % = 1.0, at which the initial static structure factor 
is given by Si(k) = S^ eq \k-,(ft,Ti) = 1 /n£h(k;(ft,Ti) (dotted lines), to a final temperature 
T = 0.2. As we can see in Figs. [5](a) and (b), in both cases the non-equilibrium stationary 
structure factor S a (k) does not exhibit any singular or non-physical feature at any wave- 
vector. Of course, for the same reason, the non-equilibrium static structure factor S(k\t ) 
will evolve smoothly from Si(k) at t — 0, to S a (k) as t —» oo, without exhibiting any hint 
of the singular behavior characteristic of the thermodynamic function \n£f(k )] _1 below the 
spinodal curve. This is illustrated in Fig. [5] by the light-gray solid lines, corresponding to 
S(k;t) evaluated at representative finite values of the waiting time t. Notice in particular 
that for the quench illustrated in Fig. [5](a) (corresponding to (ft = 0.2), a mild but noticeable 
small-/c peak of S(k;t ) develops at k max (t ) cs 1. This is to be contrasted with the quench 
along the isochore (ft = 0.4 illustrated in Fig. |5](b), where this small-h peak is imperceptible. 
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FIG. 5: Initial equilibrium static structure factor S^ eq \k; (j>, Ti ) = 1 /n£h(k] (j>, Tj) (dotted red curve) 
of the isochoric quench from the ergodic state point (</>,!)) to the state point ( </),T ) inside the 
spinodal region (panel (a) refers to the isochore (f> = 0.2 and panel (b) to 0 = 0.4). The dashed 
line represents the thermodynamic function \n£h,(k\ <j>, T)] -1 at the latter point, whose negative 
and singular features indicate the nonexistence of uniform equilibrium states with static structure 
factor S^ eq \k\^),T) = 1 /n£h(k;<f>,T). The dark (green) solid curves correspond to the (physically 
acceptable) asymptotic long-time static structure factor S a (k ) = lim^oo S(k',t) = S*{k\u a ) of this 
quench. The soft gray curves represent S(k\t) evaluated at an (arbitrary) sequence of values of 
the waiting time t. 


B. Relationship with the theory of Cahn, Hilliard, and Cook. 


One of the main general features of the non-equilibrium evolution of S(k]t) after an 
isochoric quench inside the spinodal region is precisely the development of this low-fc peak 
of S(k]t) located at the time-dependent wave-vector k max (t). This peak is associated with 
the growing length scale £(£) ~ 2'K/k max (t) of the clusters formed in the early stage of 


the process of spinodal decomposi 


Cahn, Hilliard, and Cook (CHC) 


ion. According to the classical theory developed by 


(J , and reviewed in detail by Furukawa [7j], this small- 
k peak evolves with waiting time t, moving to progressively smaller wave-vectors, so that 
Hindoo £(£) = oo. The correct interpretation, however, is that this scenario only describes 
the very early stages of the phase separation process, up to a point in which the size £(t) 
reaches mesoscopic dimensions so that additional effects (such as surface tension, convection, 
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etc.), not contained in the theory, drive the system to full phase separation. 

As it happens, the NE-SCGLE theory contains this classical CHC theory as its linear 
and small-A: limit. To see this more explicitly, let us notice that our fundamental time 
evolution equation for S(k;t ) in Eq. (l2.ip above, becomes Cook’s equation (Eq. (3.4) in Ref. 
j3]) in the linear regime and the long-wavelength limit. The linear regime is achieved in 
our theory by neglecting the non-linear dependence of the time-dependent mobility b(t) on 
S(k',t ) itself, i.e., by approximating b(t ) = 1, whereas in the long-wavelength limit (small 
fc’s) we can expand the thermodynamic function n£h(k ; 0, T) in powers of k up to quadratic 
order, n£ h (k;(j),T) & ( a/k B T) + (K/kBT)k 2 , with a = a(0, T) arid K = Ji(0, T) being wave- 
vector independent parameters. These two approximations, with M 0 = D°/k B T, convert 
Eq. (12. ip into 


dS(k; t) 


which is Eq. (3.4) of Ref. 


= 2k 2 M 0 [k B T - (a + Kk 2 )S{k ; t)] , 


(4.1) 


7]. Notice also that Cook’s approximate solution of this equation 


for shallow quenches (Eq. (3.9) of Ref. [?|) is identical to the approximation S(k;t) = 
S*(k]u(t)) ~ Si(k) + 2k 2 D 0 u(t) quoted above, provided that b{t ) = 1. 

The fact that the NE-SCGLE theory contains the CHC theory as a particular limit also 
determines the limitation of the NE-SCGLE theory to describe the full phase separation 
process. Thus, for shallow quenches, where we know that full phase separation occurs, we 
understand that the NE-SCGLE theory only describes the early stage of spinodal decom¬ 
position. Due to its non-linear nature, however, the NE-SCGLE theory complements this 
classical scenario (in which lim^^ £(t) = oo) with the prediction that Hindoo £(f) = £ a , 
i.e., with the possibility that the process of spinodal decomposition is arrested when £( t) 
reaches the length scale £ a associated with the low -k peak of the non-equilibrium asymp¬ 
totic structure factor S a (k). Of course, if the predicted value of £ a is sufficiently large for 
surface tension and convection to take over, the system will phase-separate completely. The 
possibility exists, however, that the predicted value of £ a is not that large, or to be even of 
the scale of a few particle diameters. Under such circumstances the dynamic arrest of the 
phase separation process will occur right at its early stage, freezing the current structure and 
preventing the system from phase separating completely. These conditions are illustrated 
by the quench to T = 0.2 along the isochore 0 = 0.2, whose S a (k) is represented by the solid 
line of Fig. [5](a), for which k max ~ 1/cr, and hence, £ a ~ 6cr. Thus, we need a criterion to 
locate the crossover temperature T o (0) below which the predicted dynamic arrest will in fact 
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occur, and above which full phase separation will be driven by effects not contained in our 
equations. In the rest of this section we shall search for such criterion in the temperature 
dependence of the asymptotic structure represented by S a (k). 

C. Non-equilibrium small-/c peak of S a (k] (f>,T) and diverging length scale £„(</>, T). 

Let us thus analyze in more detail the dependence of S a (k;(f),T) on the depth of the 
quench, i.e., on the final temperature T, as well as on the volume fraction qb. With this 
intention, in Fig. 13a) we present a set of results for S a (k;(j),T ) that illustrate the In¬ 
dependence of this asymptotic structural property for quenches along the isochore <p — 0.2. 
All of these quenches start with the system equilibrated at the same initial temperature 
Tj = 1.0, whose corresponding equilibrium static structure factor S'j(fc) = S^ eq \k',(j) = 
0.2, T = 1.0) is represented by the dotted curve. Each quench ends at a different final 
temperature T, and the solid lines represent the resulting S a (k](j) = 0.2,T). Thus, the 
dotted line and the solid line corresponding to S a (k](f) = 0.2, T = 0.2) are the same as the 
dotted and solid lines of Fig. [5j)a) , but now in the vertical axis we use logarithmic scale 
to visualize the strong increase in S a (k max ](j),T) as the final temperature T approaches the 
spinodal temperature T s ((j) = 0.2) = 0.76 from below. This is clearly the most visible trend 
exhibited by this set of curves, together with the fact that the position k max ((p,T) of this 
growing non-equilibrium peak of S a (k,(j),T) moves monotonically to the left, so that the 
corresponding length scale £ a (0, T) = 2n/k max ((j),T) increases with increasing T. 

In fact, to further visualize these trends, in Fig. [6](b) we plot as a solid line the height 
Sa(kmax',<f>,T ) of the low-A; peak of S a (k;(f),T) as a function of the final temperature T 
along the isochore (j) = 0.2. This curve clearly illustrates that the behavior of S a (k max ; </>, T ) 
is distinctly different in the regime where the final temperature falls below the glass-glass 
transition temperature T c (</> = 0.2) = 0.22, and in the regime where T falls in the interval 
T c < T < T s . In the former, for the lowest temperatures S a (k\ qb, T) starts almost con¬ 
stant, with a value similar to Si(k = 0) (recall that S a (k max ) ~ Si(k max ) + 2kf nax D 0 u a ~ 
Si(k max ) ~ S)(0) ~ 1), but as T approaches T c from below it increases sharply to a value 
S a {k\(j) = 0.2, T~) = 55. In the interval T c < T < T s , on the other hand, S a (k max ) starts 
from the value S a (k](j) = 0.2, T c + ) = 55 (i.e., S a (k](f),T ) is continuous at T = T c ), and in¬ 
creases monotonically up to the spinodal temperature T s , at which it diverges as ( T s — T) -/i , 
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FIG. 6: (a) Non-equilibrium stationary structure factor S a (k;4 > 1 T) as a function of wave-vector 
k for a sequence of instantaneous isochoric cooling processes along the isochore <j> = 0.2 from 
the initial temperature T) = 1.0 to the indicated final temperature T. The figure illustrates the In¬ 
dependence of the low-h peak {k ~ 1/cr; also zoomed in the inset) and of the main peak ( k ~ 2iv/a) 
of S a {k](j),T). (b) Height of the low -k peak (dark solid line) and of the main peak (dark dashed 
line) of S a (k; 4>,T), as a function of the final temperature T along the isochore 4> = 0.2. The soft 
gray solid line represents the diverging behavior (when T —>• T s from above) of the k = 0 value of 
the equilibrium static structure factor, S^ eq \k = 0; 0, T). 


with fi ps 1.75. This divergence of S a {k max 'i <f>,T) is the non-equilibrium counterpart of the 
divergence of S^ eq \k = 0; 0, T) = ^(p/ks^/dn]^ 1 when T approaches T s from above. For 
completeness, although it has no direct relevance regarding the possibility of dynamic arrest, 
the divergence of S^ eq \k = 0; 0, T) is also shown in Fig. [6](b) (soft gray solid curve to the 
right of T s ). 

Let us also notice in Fig. [61(a) that in the regime T c < T < T s , while S a {k max ) increases 
and diverges as T —> T~, its position fc maa; (0, T) moves to the left and vanishes at T — T s , so 
that the corresponding length scale £ a (0, T) = 2'jr/k rnax ((f>,T) also diverges as T approaches 
T s from below. This divergence is the non-equilibrium counterpart of the divergence of the 
equilibrium correlation length £ ei? (</>, T) as T —)■ T s + , predicted in 1914 by Ornstein and 


Zernike 


36] to explain the phenomenon of critical opalescence along the critical isochore. In 


Fig. [7] we present the temperature dependence of £„(</>, T ) in the interval T c < T < T s along 
three illustrative isochores. We notice that in this interval £ a (</», T) exhibits two sub-regimes, 
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FIG. 7: Characteristic length £ a (4>, T) = 2i T/k max ((j),T) associated with the low-A; peak of 
S a {kmax j 4>i T), as a function of the final temperature T of the quench, for the isochores (a) = 0.05, 

(b) <fi = 0.25, and (c) <fi = 0.30. The dots represent £ a (</>, T) determined directly from S a (k max ; <j), T ) 
at a discrete set of temperatures, while the composed (solid-dashed) curve is the fit of these numer¬ 
ical data obtained with the parametrization expression in Eq. (14.21) . The vertical lines represent 
the temperatures T c {<f>) < Tq{( f)) < T s (</>), with Tq{<$) indicating the crossover of the T-dependence 
of £ a (</>, T) from exponential (solid magenta portion) to diverging power-law (dashed blue portion). 


characterized by a different dependence on the final temperature. In the first sub-regime 
£ a (0, T) increases exponentially with T, whereas in the second, which we shall refer to as 
the “Ornstein-Zernike domain”, £ a (</>,T) increases with T as an inverse power of the reduced 
temperature T* = (T s — T)/T s , diverging at the spinodal line with an exponent close to the 
classical (Ornstein-Zernike) exponent that describes the critical divergence of the equilibrium 
correlation length 36]. 

In fact, we have noticed that the results for £ a (0, T) can be conveniently parametrized 
with an empiric expression that matches an exponential function below T 0 with a power law 
above T 0 , in a manner that both, the function and its derivative, are continuous at T 0 . Such 
an expression reads 


UT) 


UT c + y 


T-T c \ 

Ts-T 0 ) ^ 


any 


Tq-T c \ 

Ts-Tq) 


T s -T 


— V 


T s -T 0 


for T c < T < T 0 
for T 0 < T < T s . 


(4.2) 


Since we already have the theoretical prediction for the spinodal and the glass-glass transition 
lines T s (4>) and T c (0), and we can determine the value £ a (T c + ) of £ a (T) immediately above T c , 
this functional form allows us to determine the exponent v and the crossover temperature 
Tq(4>) by the best fit of the numerical data of £ a (T). In fact, we actually fixed the exponent 
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v to its classical OZ value v = 0.5, so that the only parameter to be determined by this fit is 
the crossover temperature To(<f>). As observed in Fig. [TJ, the quality of this pararnetrization 
of the numerical results for £ a (0, T) along the three illustrative isochores is excellent. 


D. The threshold temperature T o (0) and the gelation line. 


The procedure just described, performed at a sequence of isochores, defines a curve T o (0) 
in the dynamic arrest diagram, which divides the interval T c (0) < T < T s (0) in two well- 
defined sub-regimes: the subinterval T c (</>) < T < To (</>), where £ a (T) increases exponentially 
with T, and the subinterval T o (0) < T < T s (0), where £ a (T) enters its diverging Ornstein- 
Zernike regime. Although we cannot jump immediately to the conclusion that these two 
subintervals are to be identified, respectively, with the gel and the full spinodal decomposi¬ 
tion regions of state space, these findings provide important elements to consider in building 
a sound physical interpretation of the dynamic arrest diagram in Fig. [H The expectation, 
however, is that indeed this crossover curve will eventually be identified with the gelation 
line. 

Further elements supporting the notion of the existence of this crossover are also revealed 
by the behavior of other properties. In fact, going back to Fig. [3l one can immediately 
appreciate that at least along the isochore (j) = 0.2, also the asymptotic square localization 
length 7*(T) seems to exhibit a similar behavior as £ a (T). To test this idea, let us propose 
a pararnetrization of its T-dependence similar to that employed for £ a (T), namely, 
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7 :(T c )e° 
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Tq~t c \ 

Ts-T 0 J 



T s -T 0 


for T c < T < T 0 
for T 0 < T < T s . 


(4.3) 


We found that a reasonable representation of the numerical results for 7 *{T) in Fig. [3] could 
be obtained with the exponent a fixed as a = 1.5, and adjusting T 0 to find the best overall 
fit, thus leading to the value T 0 (</> = 0.2) = 0.52. We tested the same procedure for other 
isochores, also fixing the exponent a as a = 1.5, and determined T o (0) from the best overall 
fit. Fig. H] illustrates the accuracy of this fit at other three illustrative isochores. This 
procedure provides a second determination of the location of the crossover curve To (</>). 

Let us finally mention that not only the spatially-meaningful length scales £ a (0, T) and 
1^/7* (0, T) exhibit this dual temperature dependence. Looking back again at Fig. [31 we 
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FIG. 8 : Dynamic order parameter 7 *(</>, T) as a function of the final temperature T, and its 
parametrization according to Eq. (14.31) at the isochores (a) <^> = 0.05, (b) cf) = 0.25, and (c) 
0 = 0.30. As in Fig. [TJ the vertical lines represent the temperatures T c (</>) < To(^) < T s (0), with 
To( 0 ) indicating the crossover of the T-dependence of 7 a ( 0 , T) from exponential (solid magenta 
portion) to diverging power-law (dashed blue portion) 


see that the behavior of the “material time” u a ((J),T) also suggests a similar pattern. Al¬ 
though we do not provide further details, let us mention that we have carried out the same 
parametrization exercise, with a similar composed functional form, leading to similar accu¬ 
racy of the overall fits. Although in this case the power-law exponent does depend slightly 
on volume fraction, this parametrization with the exponent at the value of 2, does provide 
still a third route to determine the crossover temperature To (</>). Of course, there is no rea¬ 
son to expect full quantitative consistency between these three manners to determine T o (0). 
Nevertheless, they coincide in the general notion that the region bounded by the spinodal 
curve T s (0) from above and by the glass-glass transition line T s (0) from below, is divided 
in two physically distinct domains. 

E. High-density “ordinary” liquid-glass transition. 

Let us also notice from Fig. [(3(a) that besides the development of this low -k spinodal- 
decomposition peak, the asymptotic structure S a (k‘, </>, T ) also exhibits its main rneak at k & 
27r/cr, associated with the short-range nearest-neighbor arrangement among the particles. 
This main peak also varies with the final temperature T, but in a far less dramatic manner. 
Thus, we see that the height of the main peak, plotted as the dashed line of Fig. [6](b), 
increases only slightly as the final temperature increases from 0 to T c , and then decreases 
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monotonically in the interval T c < T <T S . This behavior has some form of continuation in 
the equilibrium regime (above the spinodal temperature) in the height of the main peak of 
the equilibrium static structure factor S^ eq \k] <j>,T) (soft gray solid line for T > T s ). Thus, 
if the importance of these two competing modes is measured by the height of the respective 
peak of S a (k; <j>,T), then we conclude that for shallow quenches inside the spinodal region, 
but at volume fractions below (pb, the asymptotic structure is dramatically dominated by 
the small- k peak, which diverges as we approach the spinodal line from below. Along the 
same isochores, however, deeper quenches (below the glass-glass transition temperature T c ) 
lead to an opposite scenario, in which the asymptotic structure is dominated by the main 
peak of S a (k] (p,T). With the aim of understanding the nature of this glass-glass transition, 
let us first discuss its relationship with the “ordinary” liquid-glass transition predicted to 
occur at T c ((p), but for volume fractions above the bifurcation point, (p > (pb- 

Thus, Fig. |9] presents similar results as Fig. [6j but now at the isochore (p = 0.4, represen¬ 
tative of volume fractions well above the volume fraction (pb = 0.33 of the bifurcation point. 
Clearly, the scenario is now completely different from the previous one. As we can see in 
Fig. [9](a) , the main difference is that now the small-fc peak is virtually absent. In reality 
this peak is still there, but it is imperceptible in the scale of the main figure, as revealed by 
the rather extreme zoom in the inset. Thus, in the virtual absence of this spinodal peak, 
the most visible feature of S a {k',<p,T ) is its main peak located at ka ~ 2n. As illustrated 
by the set of plots of S a (k] <p , T) in Fig. |9](a), the height of this main peak exhibits a rather 
mild variation as a function of the final temperature T. This is also illustrated by the dark 
blue dashed line of Fig. |9](b). 

Let us notice that along the isochores with volume fraction (p < (pb, the dominating 
feature is the existence of the interval T c < T <T S , where the spinodal decomposition peak 
determines the nature of the process of dynamic arrest. This interval naturally defines the 
other two complementary regions of state space: the high-temperature equilibrium domain 
above the spinodal curve, and the low-temperature glass region below T c (<p). With this 
scenario in mind, it is natural to question if at isochores with volume fraction above (pb , where 
the order T c < T s has been inverted to T s < T c , we will again find three physically distinct 
temperature regimes: the high-temperature equilibrium domain above T c (0), the window 
T s < T < T c , and the low-temperature domain below the spinodal curve T, (</>). In reality, the 
NE-SCGLE answer to this question was already given by the behavior of y a in Fig. [21 which 



FIG. 9: (a) Non-equilibrium stationary structure factor S a (k;4 > 1 T) as a function of wave-vector 
k for a sequence of instantaneous isochoric cooling processes along the isochore <j> = 0.4 from 
the initial temperature T) = 1.0 to the indicated final temperature T. The figure illustrates the In¬ 
dependence of the low-A: peak {k ~ 1/cr; also zoomed in the inset) and of the main peak ( k ~ 2iv/a) 
of S a {k](j),T). (b) Height of the low -k peak (dark solid line) and of the main peak (dark dashed 
line) of S a (k; 4>,T), as a function of the final temperature T along the isochore 4> = 0.4. The soft 
gray solid line represents the diverging behavior (when T —>• T s from above) of the k = 0 value of 
the equilibrium static structure factor, S^ eq \k = 0; 0, T). 


only identifies two physically distinct regions, namely, the high-temperature equilibrium 
liquid domain, T > T c (</>), and the region of dynamically arrested states, 0 < T < T c , 
without the spinodal temperature T s (</>) having any relevance. 

In an attempt to better understand this scenario, let us notice that along the isochores 
with volume fraction larger than (/>&, we have that the equilibrium static structure factor 
S^ eq \k] (j),T) exists in the interval T > T s (0) (S^ eq \k = O;0, T) is the soft gray solid line 
in Fig. [9](b)), and the dynamically-arrested static structure factor S a (k](p,T) exists in the 
interval T < T c (0) (see the dark (blue) solid and dashed lines in Fig. El(b)). This means 
that both stationary solutions of the time-evolution equation for S(k]t), Eq. (12.Ill , exis; 


simultaneously in the interval T s (</>) < T < T c ((j)). However, it is not difficult to see 34] 


that under these circumstances, S a (k] (j >, T) is the only possible long-time asymptotic limit of 
S(k;t). Thus, although S^ eq \k]4>,T) exists and is well-defined in the interval T s (cj)) <T< 
T c (</>), it does not have any influence in the structural relaxation of the system for quenches 
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in that interval. This is why the spinodal divergence of S^ eq \k = 0; (f),T) bears no relevance 
on the long-time arrested structure of the system, and hence, the spinodal temperature T s 
ceases to have any dynamical significance in the regime (j) > fa. 

This means that if we revisit the predicted dynamic arrest diagram in Fig. [4] we might 
as well remove completely the dashed portion of the spinodal curve to indicate its dynamic 
irrelevance and to make clear the predicted continuity of the glass phase throughout the re¬ 
gion of state space below the curve T c (0). In the high-packing/high-temperature regime this 
glass phase corresponds to the “hard-sphere” or “repulsive” high-packing glasses that result 
from compressing hard spheres to (p^ s ~ 0.58, whereas in the low-packing/low-temperature 
regime it corresponds to “attractive” glasses, in which the intense (compared with at¬ 
tractive interactions allow the formation of macroscopically homogeneous spanning clusters 
with large microscopic heterogeneities involving high local packing coexisting with voids of 
similar dimensions, to allow mean packing fractions much lower than (p^ s & 0.58. This 
predicted scenario is fully consistent, at least qualitatively, with the existence of the glass 
region of the experimental scenario reported in Refs. jl3, 14], and referred to there as region 
III. 


V. PREDICTED SCENARIO AND AGREEMENT WITH EXPERIMENT. 

In this section we shall actually revisit and update the dynamic arrest diagram of Fig. |4j 
The resulting updated diagram will constitute the core of the scenario predicted by the NE- 
SCGLE theory, whose connection with the reported experimental scenario on the arrested 
spinodal decomposition of lysozyme solutions will also be discussed. Let us thus incorporate 
the two most relevant findings of the previous section, in the dynamic arrest diagram of Fig. 
[4j The first is concerned with the dynamic irrelevance of the dashed portion of the spinodal 
curve. This leads to the predicted continuity of the glass phase throughout the region of 
state space below (and to the right of) the curve T c (</>). 

The second refers to the existence of a crossover temperature Tfafa) at which several 
properties change their temperature dependence from an exponential to a power-law that 
diverges at the spinodal line T s (fa). In Fig. [TO] we present the updated dynamic arrest 
diagram. There we have indicated the location of the crossover line T 0 (</>) obtained in the 
three manners described above (i.e., based on the T-dependence of £ a (0, T), of y*(0, T ), and 
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FIG. 10: Dynamic arrest diagram with the crossover line To(</>) determined by the parametrization 
of the transition from exponential to power law of: the typical size £ a (<j), T) of the heterogeneities 
(T°), the square localization length 7 a ((j),T) (T°), and the the material time n a ((/>,T) (T°). Each 
of these three curves provide an approximate location of the crossover from full phase separation 
(region I) to gel formation (region II). Region III corresponds to the formation attractive glasses. 

of u a ((j),T)). We do not observe, but we did not expect, any quantitative agreement among 
these three results for T o (0). However, it is clear that within their quantitative differences, 
the three of them suggest a common qualitative and semi quantitative scenario. In fact, 
their spreading may serve to remind us that we are talking about a crossover between two 
limiting conditions, and not about a sharp transition between two well-defined states (like 
in the liquid-glass transition line). 

The referred two limiting conditions are, for shallow quenches (T < T s ), the ineffec¬ 
tiveness of the dynamic arrest conditions to prevent full gas-liquid separation, whereas for 
deeper quenches, slightly above the glass-glass transition (T > T c (</>)), the effectiveness 
of the condition of dynamic arrest to interrupt the process of spinodal decomposition at 
an early stage of its development. The result of this interruption cannot be other than 
the formation of glassy states whose asymptotic structure is strongly determined by the 
amplified heterogeneities represented by the frozen small-fc peak of S a (k). These highly het¬ 
erogeneous structures contrast with the much more homogeneous structure of the attractive 
glasses formed at deeper quenches, below the glass-glass transition (T < T c (<f>)), in which 
the relevant spatial heterogeneities are determined by the main peak of S a (k), and not by 


31 







its insignificant small-A; peak. Besides the difference in the scale and degree of the spatial 
structural heterogeneity of the two glass phases, we can highlight another relevant difference, 
this time in the localization length a/ 7 *(</>, T ). From the results of the previous sections we 
can see that ^ 7 *(</>, T) has the Lindemann-like value ^ 7 *(</>, T) ps 0.15 for the glass phase 
right below T c (jp). In contrast, for the glass phase right above T c (fa), the corresponding 
value of -^/ 7 *( 0 , T) is about ten times larger, thus indicating a sudden increase in the degree 
of heterogeneity upon crossing the glass-glass transition line T c (</>). It is, of course, very 
difficult not to associate the fluffier glass phase above T c (0) with gel states, and the (rather 
loosely defined) crossover line T 0 (</>) with the corresponding gel line. 

In fact, if we now compare the topolo gy o f our updated theoretical diagram with its 
experimental counterpart in Fig. 4 of Ref. [13] (or Fig. 1(a) of Ref. jbij]), we immediately 
recognize that the qualitative mapping of one onto the other is impressive. To enhance this 
similarity with the experimental diagram, here we have also labeled, as in Refs. Q, Q, the 
three kinetically distinct regions of the theoretical diagram as I, II and III, since the inter¬ 
pretation of each of these theoretical regions, as discussed in this section, are conceptually 
consistent with the description in Refs. 13J, [l4| of the corresponding regions experimentally 
determined in the lysozyme system. The most relevant feature shared by the predicted and 
the experimental scenarios, is the fact that the “ordinary” liquid-glass transition penetrates 
the unstable region as a novel glass-glass transition. For (j) > fa, this glass phase borders 
with the equilibrium liquid phase along the liquid-glass transition line T c [fa), so that we 
expect that the predicted phenomenology of the corresponding liquid-glass transition will 
)e essentially identical to that of the purely repulsive soft-sphere system studied in Ref. 


34]. Although there is still much to be understood about this “ordinary” glass transition, 


from the study in Ref. 


34j we have learned that as soon as we enrich the scenario provided 


by the dynamic arrest diagram (obtained solely from long-time asymptotic properties like 
7 *(0, T)) with the kinetic dependence on waiting time, a number of important effects can 
be predicted, such as the aging of the structure and the dynamics, and the experimental 
impossibility of observing the sharp and discontinuous liquid-glass transition implied by the 
dynamic arrest diagram, which can only be observed as a blurred and continuous transition 
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Let us finally mention other two additional points of comparison between our theoretical 
scenario and the experimental results of Schurtenberger and collaborators [l3, lj|. This 
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FIG. 11: Comparison between the theoretical scenario and the experimental results shown in 


Figs. 3(b) and 4(c) of Ref. 


. 3 


for the characteristic length £ o (0, T) = 2-k/ k max {4>^ T). In the 


main panel of (a) we show the theoretical results for £ a as a function of T* = ( T s — T)/T s at 
4> = 0.15 for the HSAY model with z = 2, parametrized according to Eq. (14.21) . In the inset 
the dots are the experimental data in Fig. 4(c) of Ref. |l4] for the fast quench of the lysozyme 
solution at the same isochore, and the composed solid-dashed curve is the fit of these data with 
the parametrization expression in Eq. (|4.2p . yielding T 0 = 12.6°C (vertical solid line). Fig. (b) 
illustrates the experimental observation reported in Fig. 3(b) of Ref. fill ] that for fixed final 
temperature £ a (0, T) seems to be independent of volume fraction (inset). This feature is also 
qualitatively predicted by our theory, as illustrated by the results for £ a (0, T) along the isotherms 
T = 0.20,0.30,0.45 plotted in the main panel (circles, squares and diamonds, respectively). The 


experimental data in the insets were read from Figs. 3(b) and 4(c) of Ref. 
of IOP Publishing. 


. [l4] wi 


with the permission 


group measured the temperature and density dependence of the heterogeneities length scale 
£ a (0, T) of the arrested structure of the gels. They reported that for fixed volume fraction 
(0 = 0.15), £ a (0, T) could be adjusted by £ a (0, T) oc [T s — T]~ v with v = 0.55 ± 0.1 jhij ]. 
and that £ a (0, T), measured at fixed temperature, seemed to depend only weakly on volume 
fraction. Regarding the temperature dependence of £ a (0, T), in the main panel of Fig. 
m a), we present the theoretical results for £ a (0, T) in the HSAY model with z = 2 along 
the isochore 0 = 0.15 parametrized, as in Fig. [Tj. according to Eq. (14.21) . but plotted here 
as a function of the reduced temperature distance T* = (T s — T)/T s to the spinodal curve. 
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As in Fig. [7] the vertical (dotted and solid) lines represent, respectively, the glass-glass 
transition temperature T c (0) and the temperature Tq( 0) indicating the crossover of the In¬ 
dependence of £ a (0, T) from exponential (solid magenta portion) to diverging power-law 
(dashed blue portion). The prediction is that gels should be expected to form in the interval 
T c < T < T 0 and to start to phase separate for temperatures above T 0 . In the inset of Fig. 
nra) we reproduce the log-log plot of the experimental data reported in Fig. 4(c) of Ref. 
14] for the fast quench of the lysozyme solution at the isochore (p = 0.15. The composed 
(solid-dashed) curve is the parametrization expression in Eq. (14.2ft . in which we fixed the 


exponent v at its Ornstein-Zernike value v = 0.5 and employed the reported experimental 
value of the parameters T s (= 18°C), T c (= —4.5°C') and £ a (0, T c ) (= 1.05/im), and varied 
T 0 to obtain the best fit of these experimental data, yielding T 0 = 12.6°C, a value 16 % 
below the experimental gel (or tie-line) temperature of 15°C. 

The experimental observation that for fixed final temperature £ a (0, T) seems to be in¬ 
dependent of volume fraction is illustrated in the inset of Fig. fTTT b). which plots the ex¬ 


perimental data of Fig. 3(b) of Ref. [jj]]. In the main panel of our figure fTTT b) we plot 
the theoretical £ a (0, T) as a function of <p at fixed temperature, in the same volume fraction 
window as the experimental report. We observe that indeed, in this particular window, the 
theoretical prediction for £ a (0, T) along the three isotherms presented can be said to depend 
rather weakly on volume fraction. We must mention, however, that in a more extended 
window that includes lower densities, the spinodal divergence of £ a (</>, T) will lead to a more 
complex scenario than the simple independence of £ a (0, T) on volume fraction, as already 
insinuated by the curve corresponding to T = 0.45. Let us also stress that the comparison in 
Fig. fill between our theoretical predictions and the experimental results of Ref. [14] is only 
meant to illustrate the qualitative consistency between our theoretical predictions for a spe¬ 
cific model system and the results of an experimental study on the real system conceptually 
closest to our idealized model. Although the difference in the regime theoretically studied 
here (relatively long-ranged attractions, z=2) and the experimental system (much shorter- 
ranged attractions) prevents a more detailed direct comparison, it is clear that our theory 
is able to predict and describe the relevant measurable properties, and that the predictions 
are consistent with the available experimental reports. 
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VI. SUMMARY. 


In summary, in this paper we have introduced the Non-equilibrium Self-consistent Gener¬ 
alized Langevin Equation theory of irreversible relaxation in liquids, as a general theoretical 
framework for the description of the non-equilibrium processes involved in the spinodal de¬ 
composition of suddenly and deeply quenched simple fluids. The main result of the paper is 
the dynamic arrest diagram in Fig. [TO] of the previous section, which is the non-equilibrium 
analog of the equilibrium phase diagram in Fig. [0 This dynamic arrest diagram com¬ 
plements the equilibrium scenario in Fig. [Qwith the kinetic perspective of the NE-SCGLE 
theory, which predicts the dynamic arrest of a simple model liquid of attractive hard spheres 
suddenly and deeply quenched to its unstable region. Because of the non-linear and non¬ 
equilibrium nature of the phenomenon described, the initial interpretation of the specific 
predictions turned out to be rather subtle. However, the detailed analysis of several proper¬ 
ties associated with the long-time asymptotic structure factor S a (k;(j),T ) provided a rather 
coherent picture, summarized by the resulting dynamic arrest diagram of Fig. [TO] 

The first main feature of the predicted dynamic arrest diagram of attractive hard spheres 
is the continuity of Region III, corresponding to the glass phase below and to the right of 
the curve T c (0), which contains the high-packing/high-temperature “repulsive” glasses and 
the low-packing/low-temperature “attractive” glasses. For 0 > 0&, this glass phase borders 
with the equilibrium liquid phase along the liquid-glass transition line T c (0). For 0 < 0& 
this liquid-glass transition penetrates the unstable region of the gas-liquid coexistence as a 
glass-gel transition, with the gel phase constituting region II of the diagram. Finally, we 
associated the formation of gels with the early arrest of the growing spatial heterogeneities, 
whose size £ a (<fr,T) = 2n/k max is determined by the position k max of the spinodal (small-/;;) 
peak of S a (k](f>,T). The temperature dependence of this and other properties revealed the 
existence of two regimes, separated by a crossover temperature T o (0). This crossover line 
was conjectured to separate the region of gel formation from the region where the dynamic 
arrest is not able to prevent full phase separation. 

This theoretical scenario emerged only form the analysis of long-time asymptotic struc¬ 
tural properties, such as S a (k]<j>,T). The NE-SCGLE theory, however, provides abundant 
additional information, including the detailed non-equilibrium evolution of the static and 
dynamic properties of the system at any instant after the quench. This kinetic and dynamic 
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information leads to a richer and more detailed scenario whose description, however, is out 
of the scope of the present paper. Nevertheless, even without this important information, 
we related the available theoretical scenario with the results of the study of dynamically 
arrested spinodal decomposition in the experimental model system conceptually closest to 
the theoretical model analyzed here. As we demonstrated in the previous section, the most 
relevant predicted features are in satisfactory qualitative agreement with the experimental 
observations. 

Of course, rather than closing the theoretical discussion of arrested spinodal decomposi¬ 
tion, the results presented in this paper only open a number of question to further discussion. 
For example, we expect to expand this discussion to include the kinetic aspects not discussed 
here, as well as the dependence of the resulting dynamic arrest scenario on the main fea¬ 
tures of the interaction potential, such as its spatial range and its specific form (attractive 
Yukawa vs. square well, hard-sphere vs. soft sphere, Lennard-Jones, etc.). Similarly, we 
have not exhausted all the pertinent tests and comparisons of the theoretical predictions 


with available re 


evant simulations, such as the thorough simulation study of Testard, Kob 


and Berthicr 19]. In future communications, however, we shall address these and other 


issues that we could not properly treat in the present paper. 
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Appendix. The NE-SCGLE theory: main approximations. 

Let us have in mind a simple monocomponent liquid formed by N Brownian particles 
in a volume V, which diffuse with a self-diffusion coefficient D° between collisions while 
interacting through a pair potential u(r). The NE-SCGLE theory was originally derived 
from a non-equilibrium extension of Onsager’s theory of thermal fluctuations [32]. This 


36 







abstract theory provides the general structure of the time-evolution equations for the mean 
value and the covariance of the thermal fluctuations of the macroscopic variables of a non¬ 
equilibrium system. Applied to the local number concentration of particles in our Brownian 
liquid, these two equations become 

( --J— = D ° V • t) W/4 r ; «(*)] (!) 

and 

=D°V-w(r,t) 6(r,t)V J dr 2 £[r, r 2 ; n(t)]a(r 2 , r'; t) 

+ D°V' -n(r',t) b(r',t)V I dv 2 £[r\ r 2 ; n(f)]cr(r 2 , r; f) ( 2 ) 

- 2D°V • n(r, t) 6(r, t)V5(r - r'). 


which describe the non-equilibrium evolution of the mean value n(r,t) and the covariance 


cr(r, r'; t) = 5n{ r, t)Sn( r', t) of the fluctuations 5n( r, t) = n(r, t) —n( r, t) of the local concen¬ 
tration prohle n(r, t) of a colloidal liquid. 

There are two fundamental elements entering in these equations. The first involves the 
relevant thermodynamic information, and is the electrochemical potential /x[r;n(i)] of one 
colloidal particle at position r, from which the thermodynamic stability function £[r, r'; n\ = 
[<5/3p,[r; n]/Sn( r')] derives jl, 3] [ksT = 1 / fd being the thermal energy per particle). The other 
is the time-dependent local mobility function b(r,t), which in reality is also a functional of 
n(r,t) and cr(r, r';f), thus introducing a non-linear character to these equations. In what 
follows we provide a more detailed definition of these objects. 


1. Thermodynamic elements. 

The electrochemical potential /3/i[r; n{t)} is an ordinary function of the particle’s position 
r and a functional of the local concentration prohle n(r,t). Although not explicitly stated, 
/i[r;n(£)] is also a functional of the temperature prohle T(r, t). Here, however, we shall 
only have in mind uniform temperature prohles, T(r, t) = T(i), which will play the role 
of the external control parameter describing the imposed thermal treatment of the system. 
The functional dependence of fdfi[r ; n] on the arbitrary prohle n( r) embodies the chemical 
equation of state of the system, and is a fundamental thermodynamic input of the present 
theory. The functional derivative £[r, rn] = [A/5/x[r; n\/5n{r')\ is in general an ordinary 
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da(k] r, t) 
dt 


function of both, r and r', and a functional of n(r,t). Evaluated in the absence of external 
fields and at an uniform profile n(r) = no, this two-point thermodynamic property only 
depends on the distance |r — r'| and its functional dependence on the profile n(r) becomes 
an ordinary dependence of the scalar argument no, i.e., £[r, r';n] = £h (|r — r'|; no). The so- 
called direct correlation function c(r; no) is defined in terms of this thermodynamic function 
by the relationship £ h (r',n 0 ) = <5(r)/n 0 — c(r;n 0 ). 

In order to simplify the analysis of the equations above, at this point we introduce 
the “local-homogeneity” approximation, in which we approximate this two-point thermo¬ 
dynamic property by £[r, r';n(f)] = £ h (|r — r'|;n(r, t)), i.e., by the thermodynamic matrix 
evaluated at an uniform concentration profile, whose constant value is the local and instan¬ 
taneous concentration n( r, t). Under these conditions we expect that the covariance can also 
be approximated as cr(r, r'; t) ~ a( |r — r'|; r, t), so that Eq. (J2J) can also be written as 

= — 2k 2 D°b(r, t)n( r, t)£h(k ; n( r, t))a(k ; r, t) 

(3) 

+ 2k 2 D°b(r, t)n( r, t), 

where cr(fc;r, t) = f d 3 /ce lk ' x cr(|x|; r, t) and £ h (k;n(r,t)) = (2 tt)~ 3 f d 3 ke~' lk ' x £ h (\x\;n(r,t)). 

To actually solve these equations, we must indicate the specific nature of the system, i.e., 
the specific inter-particle interactions, embodied in the pair interaction potential u(r), and 
the applied external force F(r) = — acting on a particle at position r, whose effects 

are introduced through the chemical equation of state. This is written in general as 37], 

/3/x[r; n] = /?//(/?) + lnn(r) + /70(r) - c[r; n\. (4) 

The first two terms on the right side are the ideal gas contribution to the chemical po¬ 
tential, whereas the term — c[r; n] contains the deviations from ideal behavior due to inter¬ 
particle interactions. Thus, the functional dependence of c[r; n\ on the concentration pro¬ 
file n(r) must be specified. One possibility is to propose a theoretical approximation for 
this dependence, which in the language of density functional theory [37] is actually equiv¬ 
alent to proposing an approximate free energy functional. We quote, for example, the 
simplest such approximation, referred to as the Debye-Hiickel or random phase approxima¬ 
tion, in which c[r; n\ is written as c^ RPA ^[r; n\ = —f3 J d 3 r'u(\r — r'DnfV). An approximate 
c[r; n] leads to a corresponding approximate thermodynamic matrix £[r, r'; n\. For example, 
£( RPA ' > [r,r';n\ = 5(r — r')/n(r) + /3u( |r — r'|). Other more elaborate approximations are 


available, however, in the literature of the equilibrium theory of inhomogeneous fluids 


37]. 
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We have assumed so far that the external fields are static and represented by ^(r). There 
is, however, no fundamental reason why we have to restrict ourselves to these conditions. 
In fact, the general equations of the NE-SCGLE theory above can be used, within the 
range of validity of the underlying assumptions, to describe the response of the system to 
prescribed time-dependent external fields -0(r, t) or to programmed thermal or mechanical 
constraints described by a time-dependent temperature T(t) of the reservoir and/or an 
imposed time-dependent total volume V(t) of the system. Here we shall assume, however, 
that such time-dependent fields and constraints could have been used to drive the system 
to a prescribed initial state with mean concentration profile n°(r) and covariance cr°(/c;r), 
but that afterward the external fields are set constant, = ^(r). 


2. Absence of external fields and spatial heterogeneities. 


The NE-SCGLE Eqs. (JT[) and (j3]) are supposed to predict how the system relaxes to its 
final equilibrium state whose mean profile and covariance are n eq ( r) and a eq (k ; r). Describing 
this response at the level of the mean local concentration profile n( r, Dis precisely the aim 
of the recently-developed dynamic density functional theory (DDFT) 38|, |39|, whose central 
equation is recovered from our theory by setting b(r,t) = 1 in Eq. (JT]). DDFT has been 
applied to a variety of systems, including the description of the irreversible sedimentation of 
real and simulated colloidal suspensions 40]. Similarly, the time-evolution equation for the 
covariance, within the additional small-wave-vector approximation Sh(k) ~ So + ^ 2 ^ 2 + S 4 k 4 , 
can be recognized as the kinetic equation that governs the irreversible evolution of the static 
structure factor of a suddenly quenched liquid, developed in the description of the early 
stages of spinodal decomposition (see, for example, Eq. (2.11) of Ref. [8j], in which S 4 = 0, 
or Eq. (23) of Ref. [h]]). In fact, the main motivation of the present work is to contribute 
to the development of the theory of spinodal decomposition starting from the general NE- 
SCGLE theory just introduced. 

The most essential aspect of the irreversible process of spinodal decomposition is the 
development of strong spatial heterogeneities. These originate in the thermodynamic spin¬ 
odal instability of thermal fluctuations of certain wavelengths, whose amplification generates 
strong spatial heterogeneities, leading eventually to the full separation of the system in two 
separate regions. At least in its initial regime, the main features of this irreversible phase 
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separation process must be exhibited by the time-evolution of the two observable variables 
involved in the present theory, namely, the mean value n( r, t) and the covariance a(k\ r, t), as 
recognized since the early work of Cahn and Hilliard [jj] and Cook 6]. Our NE-SCGLE equa¬ 
tions actually contain Cook’s time-evolution equation for the non-stationary static structure 
factor S(k]t), but goes beyond this classical theory by incorporating the possibility of dy¬ 
namic arrest. 

To see this, let us imagine that we subject our system to a prescribed protocol of thermal 
and/or mechanical treatment, described by the chosen temporal variation T(t) and V(t) of 
the reservoir temperature and total volume, but in the absence of external fields, ^(r, t) = 0. 
We also assume that the thermal diffusivity of the system is sufficiently large that any tem¬ 
perature gradient is dissipated almost instantly (compared with the relaxation of chemical 
potential gradients), so that the temperature field inside the system remains uniform and 
instantly equal to the imposed temperature of the reservoir, T(r, t) = T(t). 

In contrast, since the essence of spinodal decomposition is the development of strong 
spatial heterogeneities, in this case we cannot in general assume that the mass diffusivity 
of the system is sufficiently large that any gradient of chemical potential can be dissipated 
almost instantly. If this were the case, the density field n(r, t) inside the system would 
always remain uniform and instantly equal to the imposed bulk density, i.e., n( r, t) = n(t ) = 
N/V(t). Under such conditions, rather than solving Eqs. ([Tj) and (]3]) for n( r, t) and a(k ; r, t), 
we would only have to solve Eq. (j3J) for a[k\ r, t), with n{ r, t) substituted by its imposed bulk 
value n(t), which now becomes a control parameter. Although the resulting scheme might 
seem too gross for the description of the spontaneous spatial heterogeneities represented by 
the neglected deviations An(r,f) = n(r, t) — n(t ) (which are definitely non-zero), here we 
shall see that it provides a welcome mathematical and numerical simplification of Eq. (J3]) 
for cr(k]r,t) ~ a(k] t), which may now be written as 

da(k\ t ) 


dt 


= —2k~D°n(t)b(t)£h{k]n(t))a(k]t) + 2k 2 D 0 n(t) b(t). 


(5) 
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3. Determination of the t-dependent mobility b(t). 


As explained in Refs. 32] and [33], the time-dependent mobility b{t) is in reality a 
functional of a(k]t), and is determined by the generalized Einstein relation 


b(t ) = [1 + / drA(*(r; t)] 
Jo 


-i 


( 6 ) 


with the f-evolving, r-dependent friction function A£*(r;t) representing the effects of the 
interparticle interactions on the total friction coefficient of a representative tracer particle. 
The functional dependence of A£*(r;i) on <j(k]t) introduces the essential non-linearities of 
the problem, and this functional dependence can only be determined in an approximate 
manner. 

If we were dealing with an equilibrium system, we could simply adopt the correspond¬ 
ing approximation proposed by mode coupling theory or by the equilibrium version of the 
SCGLE theory. The non-equilibrium SCGLE theory extends the latter to non-equilibrium 
conditions, thus leading to an approximate expression for A£*(r;i), namely, 


A C*(r;t) = 


Dn 


elk k 2 


S(k't) 


S(k; t) 


(7) 


247r 3 n(f) J 

xF(k,T;t)F s (k,r-t ) 

in terms of the non-stationary static structure factor S(k',t ) = a(k;t)/n(t ) and of the col¬ 


lective and sel 
In Ref. 


non-equilibrium intermediate scattering functions F(/c,r;t) and Fg(fc,r;t). 


32] it was indicated how to extend the exact memory-function expressions for 


F eq (k , r) and eq Fs(k , r) to non-equilibrium conditions, thus leading to the corresponding 
exact memory-function expressions for F(k,T]t) and Fs(k,T',t). It is at the level of the 
corresponding memory functions where one has to introduce the essential approximations, 
also extending the arguments of the SCGLE theory to non-equilibrium conditions. This con¬ 
verts the exact memory function expressions for F(k, r; t ) and F$(k , r; t) into the following, 
approximate memory-function equations, written in terms of the Laplace transforms (LT) 
F(k, r; t ) and Fs(k, r; t), as 

S(k; t) 


and 


F(k,z;t) = 


F s {k,z;t) = 



k 2 D°S~ 1 (k,t) ’ 

1+A(fc) AC*(«;*) 


1 

z + 

k 2 D° 


( 8 ) 


(9) 
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30], given by 


with A (k) being a phenomenological “interpolating function” 

\(k) = l/[l + (k/k c ) 2 }, (10) 

where the phenomenological cut-off wave-vector k c depends on the system considered. In 
the particular case of an instantaneous isochoric quench, which is the subject of the present 
paper, these equations become Eqs. (I2.ip ~ fj2.6p of Sect. HU 
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